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The factor of lo-'® is put in (7) because values of R„ are expressed in square nanoteslas. If 
tile root mean square value of the dipole moment per unit area is divided by the thickness of 
the magnetized layer, we shall have obtained a root mean square value of the magnetization. 
We have summed the coefficients from = U to = 23 by using the data from 
Langel & Estes i iqSa) Ibr an oceanic crustal layer ti km thick, and obtained a magnetization of 
l.t>.5x lii-®e.m.u. m'* for the magnetization of an oceanic crustal layer. This mag- 

netization is only for harmonics between 14 and 23. If the 'rrustal' portion of lower degree 
harmonics {Langel & Estes 1982) and a'l the higher harmonics from 23 up to harmonics 
equivalent to the wavelength of the sea-floor spreading anomalies are included, a much higher 
magnetization would be obtained. Typical magnetizations of oceanic crustal rocks .appear to 
lie much too low for such magnetic potentials to be generated 'Harrison 1976. 1981). 

In these three examples (equations iti) and (7) and figure 1) we have shown that the ‘core’ 
field as generated by a certain number of low degrees of spherical harmonic does not com- 
pletely remove long-wavelength terms from the surface scalar field. It is also possible to show 
that the same is true for field components, which follows simply from the higher degree terms 
in the associated Legendre polynomials. These all contain powers of cos 0 going from powers of 
zero or one (depending on the degree and order of^he harmonic) up to higher powen. We 
therefore suggest that a better model of the core field is to be found by modelling it with 
non-central sources, such as current loops located within the core. 

Detailed analyses of this nature have been done by .\lldredge & Hurwitz ^964) and 
.Alldredge & Stearns (t969). .Alldredge & Hurwitz (1964) found that to model the non-dipx)le 
field accurately, the oflf-centred dipoles needed to be placed deep within the Earth’s core. The 
probable reason for this is that the dipoles are an attempt to model a source of finite dimensions 
at the core-mantle boundars’. 

The use of radial dipjoles allows us to model the long-wavelength portion of the field better 
than the use of spherical harmonics. Each dipole can be modelled as a series of harmonics, 
which, if the origin of the spherical coordinate system is taken where the a.xis of the dipxtle 
impinges on the surface of the Earth, are all zonal harmonics. .Alldrcdge (1980) has produced 
formulae for these harmonics. It is then possible to plot the mean square field generated by a 
dipclc as a function of the degree of zonal har^monic, for different radial distances of the 
dipole from the centre of the Earth. Figure 2 shows three such plots. One is for dipoles located 
at the surface of the core. The other two are for dipoles placed at the maximum and minimum 
distances estimated by .\lldredge & Hurwitz 1964). .Also shown is the wa\ in which the 
spherical harmonics of the Earth fall off, taken from Langel & Estes 1982 It can be seen that 
each dipole, if buried fairly deeply within tlic core, as required liy the models of .Alldredge 
Hurwitz <1904) contributes harmonics which in general fall ofl'more steeply than the spherical 
harmonic representation of the Earth's field. 


Off-centred dipoles and the field in Iceland 
Cox (1975) suggested that there was a certain pattern of non-dipole sources that explained 
some puzzling observations that he had made of excursions seen in lava flow data from Havsaii. 
These observations were that the excursiems seemed .ilwavs to bi in the direction such that 
shallower inclinations of the field vserc produced. Ciox suggested tl at this was caused b> core- 
surface sources drifting past the longiuidi of Havsaii. which vs ere pi<*ferenti,illy orieiiied and 
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Modelling the core magnetic field of the Earth 


By C. G. a. Harrison and H. M. Carle 
Rosenstiel School of Marine and Atmospheric Science, University of .Miami, 

4*100 Rickenbacker Caitsetvay, .Miami, Florida 33 149, U.S.A. 

The magnetic field generated in the core of the Earth b often represented by spherical 
harmonics of the magnetic potential. It has been found from looking at the equations 
of spherical harmonics, and from studying the values of the spherical harmonic 
coefficients derived from data from Magsat, that this b an unsatisfactory way of 
representing the core field. Harmonics of high degree are characterized by generally 
shorter wavelength expressions on the surface of the Earth, but also contain very long 
wavelength features as well. Thus if it is thought that the higher degree harmonics 
arc produced by magnetizations within the crust of the Earth, these magnetizations 
have to be capable of producing very long wavelength signab. Since it is impossible to 
produce very long wavelength signab of sufficient amplitude by using crustal 
magnetizations of reasonable intensity, the separation of core and crustal sources by 
using spherical harmonics is not ideal. We suggest that a better way is to use radial 
off-centre dipoles located within the core of the Earth. These have several advantages. 
Firstly, they can be thought of as modelling real physical current systems within the 
core of the Earth. Secondly, it can be shown that off-centred dipoles, if located deep 
within the core, arc more effective at removing long w-avelcngth signals of potential or 
field than can be achieved by using spherical harmonics. The dbadvantage b that it 
b much more difficult to compute the positions and strengths of the off-centred dipole 
fields, and much less easy to manipulate their effects (such as upward and downward 
continuation). But we believe, along with Cox and Alldredge & Hurwitz, that the 
understanding that we might obtain of the Earth's magnetic field by using physically 
reasonable models rather than mathematically convenient models is very important. 
We dbcuss some of the radial dipole models that have been proposed for the non- 
dipole portion of the Earth’s field to arrive at a model that agrees with observations 
of secular variation and excursions. 
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Introduction 

The easiest way of describing observations of the magnetic field obtained over the surface of 
a sphere is to represent the magnetic potential as a sum of spherical harmonic functions of the 

r = n V fa/r)"-' (.gll' cos + a;;* sin I cos <>),' • ' (1) 


' I lu-O 


where 1' is the magnetic potential, n is the degree and m the order of the harmonic, a is a 
reference radius usually the radius of the Earth), r b the radius at which the observation takes 
place, 6 is the colatitude, ^is the longitude, (cos 0) are the Schmidt semi-normalized 
associated Legtfndrc [xilynomials, and o"' and A"' are the Gauss coefficients. Equation (1) 
assumes that there arc no sources external to the radius^ If 1 1) b differentiated along spatial 
coordinates, it gives the magnitude of the field components along the coordinates. Thus if 
Fs, F^ and l\ are the north, east and vertically downward components of field, they can be 
(wpressed as tollows; 
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It is possible to express tlie power in each dcijree of liarmonic by using the formula of Lowes 
(iq66). He showed tiiat the mean square value oxer the surface of the sphere of the field H 
produced bv harmor'cs of a given degree n. called /?,, is given by the following expression: 



/?» = («+») 1 {(.«?)= + (AD*}. (3) 

HI 

Langel & Estes ■ i«)82 ) have produced the most accur \te set of harmonics for the Earth's 
magnetic field, up to degree 23. They used data collected bv M<igsat on a scries of quiet days 
during the lifetime of this satellite. If Ig R„ from this model s plotted against n, it turns out that 
the values between n = 2 and w = 13 fall on a steeply sloping line of negative slope, whereas 
the values for n > 13 fall on anodier line with a very small negative slope. The value for 
n = 1 demonstrates that the dipole term of the field is considerably more important than any 
other. 

Bullard {1967) correlated the maximum degree and order of a set of spherical harmonics 
with a short wavelength limit of what is seen in a two-dimensional portion of the Earth’s 
surface, and showed that spherical harmonics truncated at degree and order n will have the 
same short-wavelength cutoff as a two-dimensional Fourier series of the form 

S S fl/K, cos ( 2tt/ L) ( px 4. gy) 



for which + < nL/c, where L is the length of the side of the square area and c is the 

circumference of the Earth. It is necessary to separate fields into core and crustal portions for 
two reasons. If we wish to look at anomalous magnetizations within the crust of the Earth 
generated by remanent or induced magnetization of ferromagnetic particles below their Curie 
point, then it is necessary to remove the field generated within the core of the Earth, w hereas if 
we are interested in looking for sources of the core field, it is conversely necessary to filter out 
that component generated within the crust. Since these components should have grossly 
different wavelengths, Bullard suggested that the cutoff in spherical harmonics should be 
chosen from the standpoint of the maximum or minimum wavelength that it is desired to have 
within the spherical harmonic expansion. 

There arc two problems with this method. The first is that we do not measure potential at 
the Earth's surface, but rather either a component of the field (such as the vertical component), 
or more commonly the magnitude of the total field vector. The second problem is that higher 
harmonics of the potential have within them very long wavelengths, so that the separation 
between long and short w avelengths is not perfect. NVe shall show that if total fields are being 
used the first problem means tliat wavelengths equivalent to c/'2n arc cont.'iined within 
harmonics up to degree n, meaning that shorter wavelengths are produced. W'e shall also discuss 
the second problem at some length. 

On looking at the first problem in more detail, if the observed field is in fact the scalar total 
field, then in order to generate this from the spherical harmonic^ (s.h.) potential, it is necessary 
to obtain the three spatial derivatives as in 2', and then to square and add them. The process 
of squaring effectively halves the minimum wavelength represented by the harmonics. The 
process of taking the square root to obtain tlie scalar field also adds complic,itions to the 
process, with the result that we have sometimes used the square of the total scalar lieid to 
reduce this one complication. 
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\V> shall show that the second problem is csseniialls unsolvable il's.li. are used to describe 
the held. Rather, ue simt?est that the core field shmdd be represented by spatially distinct 
sources within the core. Thus we I'ollow C ox liqbS). .Uldred^e & Hurwitz ig 64 ) and others 
who have sui^gested that diprrles located within the core may oMer ijreater insight into the / 
mechanisms w hereby the core field is generated. We show that the use of radial dipoles buried 
within the core to represent core surface current ItMrps can lead to a better approximation to 
the longer wavelength |>ortinns of the field than does the s.h. representation of the magnetic 
potential. 

Data from Iceland and other areas with good secular variation records can l>e used to 
determine something about the nature of the source regions of the non-dipole field. The 
Icelandic data are especially important in that Iceland is at a very high latitude, and thus 
serves to check various models of excursions and reversals that have been proposed. The 
Icelandic data show several differences from data from lower latitudes. FirstK , the frequency 
of virtual geomagnetic poles (v.g.ps) of low latitude is considerably greater here than else- 
where. Secondly, the longitude of v.g.ps of low latitude from Iceland appears to be random, 
compared with longitudes of low-latitude v.g.ps from lower latitude sites. In this latter case, 
when the v.g.ps are from records of reversals, the longitude appears to be confined to lying close 
to the longitude of the site. For these and other reasons, Dodson igSo) suggested that the 
non-dipole field could be modelled by radial dipoles located within the core, in which the 
probability of finding a dipol • within a certain area of core surface increased towards the poles. 

This is another example of using non-central dipoles to model the non-dipole field. 


T.\Bt.E 1. .\MPl.ITtnES OF FOfRIER COEFFICIENTS I 
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Sep.\r.\tion of the field into core .\ND CRCSTAL COMPONENTn 
Suppose we take the spherical harmonic representation of the first three zonal harmonics of 
the field by using ,1). The result is 

r = a[(rt/r - o\* cos + 'fl/r'* I cos- fi - \ ) + a/r)* j;j! I .) cos - •< cos i 4) 

Equation 2 can now be used to determine any one component of the field by differentiation 
in the appropriate spatial direction. Of more use to us is to obtain an expression for the scalar 
field, as this is the quantity that is most commonly measured when magnetic field surveys are 
made. This requires squaring, adding and then taking the square root of the three spatial 
derivatives of 4 The process of taking the square rrKit complicates matters vshrn vve trv to 
oredict w hat a simple Fourier transform of a linear profile would look like C arle Harrison 
1982 , and so we have decided to work vuth the squared scalar field. .\n expression for this is 
given Irelovv : 

T‘ » a / 1 '• - 1 cos* 0 1 -f 1 1 * a / 1 ' .j,'®,?" c^os ii/n'* ej ;i|j 1 J.l cos^^ - u o>s -d — 

*1 '*/'■ ' .»? ■ 4 .) cos — 1 X ens -0 -4 !>) 4- a/r * c" "® 4.) cos - :{o *■ *> COs d 

1 7 •> cos — I li.') 1 1 IS *0 — cos . -i 


* j a/r 
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li is convenient to express this result in terms of cosines of integral multiples of 0 rather than 
in terms of powers of cosines of y. The resulting expression is shown l)elow. This has been 
calculated for the surface of the Earth, w here r = n: 


U 1 

(«) 

Just as an example, suppose that we imagine that the ‘core* field is represented by the first 
two zonal harmonics, and that the crustal field is represented by the third harmonic. Removal 
of the 'core' field will then produce an e.xpression consisting of the last three terms in (6). The 
‘resulting 'crustal' field will have terms in integral multiples of fi all the way from zero to six. 

If wt now take a longitudinal profile of this zonal field and perform a Fourier analysis, each 
expression in the last three terms of will give a discrete Fourier component whose wave- 
length is equal to the circumference of the Earth divided by the integral multiplier of 0. Thus 
there will be Fourier components with wavelengths all the way from c/0 to r/b. The values of 
the amplitudes of these terms are given in table I, by using values of the zonal harmonics from 
a recent spherical harmonic analysis of the field recorded by Magsat. It is immediately obvious 
from this table that the terms resulting in the interaction of the first and third harmonics are 
the largest, and have wavelengths of r/U, c/2 and r/ 4 . Other wavelengths ar^ alsojmportant. 

When it is realized that the values shown in the table are in units ofX!>>’^'f^) we can see that ^ 
amplitude of the c/o harmonic is about 24 that of the r /4 harmonic. The result of the inter- | 

action on an odd harmonic (5®) with the dominant first degree harmonic ( «!') is to produce ' 

Fourier components in which the even harmonics predominate. It would perhaps be better to 
use a spherical harmonic model in which all harmonics zonal, tesscral and sectoral) are 
included. But the problem is that to write these terms down in order to illustrate principles is 
e.xtremely tedious and lengthy. Zmuda ( 1958) h- ritten down the total field squared produced 
by all harmonics up to degree 2. The resulting expression consisted of 21 terms. 

Obviously, the above example is not particularly realistic. .\ more realistic example would 
be to take a complete spherical harmonic representation of the magnetic potential, including 
core and crustal sources, and to work out total field minus core licld for a series of profiles along j 
great circles, using as a model for the core field that field generated by tlic first n degrees of 
spherical harmonic, with n about IS. Unfortunately, such a model does not exist. It would 
require about l.rt x 10^ coefficients to achieve a resolution of 100 km, which is still ‘too large to 
show up many crustal anomalies f,-\lldredge iqSt'. However, we may get an inkling as to what 
might happen by looking at the degree 23 spherical harmonic analysis produced by Langel & 

Estes ( 1982 in which they estimated that ‘core’ sources are important up to degree IS and that 
‘crustal’ sources are important for hat monies of degree l."> or above. We sliall consider the 
core' field to be caused by harmonics up to degree IS, and shall produce Fourier analyses of ^ 

residual or anomalous fields along great circle paths. These Fourier spectra are shown in 
figure 1. 


^ i + -3 cos 2ff) + S oj S cos 0 + cos Sfl) 

+ 32 [g2)~ I • ‘lo + 1 cos 2 ^ -t- 4.1 cos 
co&-2^ -t- 2.1 cos 4 ^)-' 

4- 27 -r Tii cos 2fJ + 2.1 cos 4<^) 

+ cos 0 -f- 1 Do cos :id + 4.1 cos of)) 

+ -1-70.1 cos 26 + S90 cos 46 + 1 75 cos li^}. 
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These Fourier spectra show a plateau of power running from the zero harmonic to about 
the -3rd harmonic, after which the power decreases until round-off errors produce a spurious 
signal at about harmonic number 33. It is possible to predict what the power spectrum would 
be from a hypothetical profile generated from a model in which higher degree liarmonics were 
present. It would be expected that the power in the first 24 harmonics would be higher, 
because each additional degree of spherical harmonic used in the model to describe the crustal 
field would add power in this region. But we would also obtain much more power in harmonics 
l)etwe< n 24 and the maximum degree of the spherical harmonic model, after which the power 
would fall of!' rapidly to low levels. Figure 1 emphasizes that the removal of low-degree spherical 
harmonics, considered to Ije a ‘core’ field, does not remove all power at long wavelengths from 
surface scalar fields. 



harmonic < 

Fiovre 1 . Fourier spectra of the residual field obtained by subtracting the scalar field generated bs the first 
13 degrees of spherical harmonic from the scalar field generated bs the 23 degrees of spherical harmonic. 
The 23 degree harmonic field is that described bs Langel & Estes (igSa). ftil The profile is along the prime 
meridian with a fundamental length of the circumference of the Earth. The profile is along the equator. 


It is also possible to determine directly the magnetization necessary to produce the harmonics 
of potential thought to be caused by crustal magnetization. Chapman & Bartels 1940) have 
shown that there is a direct relation between the spherical harmonics of a vertically polarized 
thin magnetic shell and the spherical iiarmonics of its potential. Tlie vertically magnetized shell 
is characterized by a vertical dipole moment per unit area. To obtain the Gauss coefficients for 
this function, it is necessary to multiply the <j"; and A'" from (I) by the quantity 


a / 2« -i- I 
4n ' II 



.Mso, from Chapman & Bartels, the mean square v alue of P"' cos O' cos Im^i ov er the surface 
t<f the sphere is equal to 1/ 2n-t- 1 . Therefore, from the values of R„ 1 equation i3 given by 
Langel & Estes 1982 we ran derive the mean square value of the vertical dijxtle moment 
per unit area: 

( 2n -t- I ) /?„ / a , a 


f 


11 


original page « 

Or POOR QUAUTT 

tiOil 


C. G. A. HARRISON AND H. M. CARLE 



I 

1 

I 





) 



I 

I 


f\l\rr>''^ n'i(. 



The factor of 10-'® is put in ( 7 ) l)ecause values of R„ arc expressed in square nanotesias. If 
the root mean square value of the dipole moment per unit area is divided by the thickness of 
the mai^netized layer, we shall have obtained a root mean square value of the mai^netization. 
We have summed the coefficients from = U to = 23 by using the data from 
Langel & Estes iiqSa) lor an oceanic crustal layer t> km thick, and obtained a magnetization of 
1 . 0.5 X 10-®e.ni.u. il.or>.\ m'’> for the magnetization of an oceanic crustal layer. This mag- 
netization is only for harmonics between 14 and 23 . If the 'crustar portion of lower degree 
harmonics (Langel & Estes 1982) and ah the higher harmonics from 23 u{ to harmonics 
equivalent to the wavelength of the sea-floor spreading anomalies arc included, a much higher 
magnetization would be obtained. Typical magnetizations of oceanic crustal rocks appear to 
l)e much too low for such magnetic potentials to be generated Harrison 1976, t98tj. 

In these three examples (equations 18) and ( 7 ) and figure I) we have shown that the ‘core’ 
field as generated by a certain number of low degrees of spherical harmonic does not com- 
pletely remove Iong-wa\elength terms from the surface scalar field. It is also possible to show 
that the same is true for field components, which follows simply from the higher degree terms 
in the associated Legendre polynomials. These all contain powers of cos 0 going from powers of 
zero or one (depending on the degree and order of^he harmonic) up to higher powers. We 
therefore suggest that a better model of the core field is to be found by modelling it with 
non-central sources, such as current loops located within the core. 

Detailed analyses of this nature have been done by .\lldredge & Hurwitz ^964) and 
.Alldredge & Stearns (1969). .Alldredge & Hurwitz (1964) found that to model the non-dipole 
field accurately, the off-centred dipoles needed to be placed deep within the Earth's core. The 
probable reason for this is that the dipoles are an attempt to model a source of finite dimensions 
at the core-mantle boundary. 

The use of radial dipoles allows us to model the hjng-wavelength portion of the field better 
than the use of spherical harmonics. Each dipole can be modelled as a series of harmonics, 
which, if the origin of the spherical coordinate system is taken where the a,xis of the dipole 
impinges on the surface of the Earth, are all zonal harmonics. .Alldrcdge (1980) has produced 
formulae for these harmonics. It is then possible to plot the mean square field generated by a 
dipclc as a function of the degree of zonal har^monic, for different radial distances of the 
dipole from the centre of the Earth. Figure 2 shows three such plots. One is for dipoles located 
at the surface of the core. The other two are for dipoles placed at the maximum and minimum 
distances estimated by .\lldredge & Hurwitz 1964J. .Also shown is the was in which the 
spherical harmonics of the Earth fall ofl", taken from Langel & Estes 1982 . It can be seen that 
each dipole, if buried fairly deeply within the core, as required by the models of .Alldrcdge & 
Hurwitz (1964,' contributes harmonics which in general fall ofl more steeply than the spherical 
harmonic representation of the Earth s field. 


Off-centred dipoles and thf field in Iceland 

Cox (1975) suggested that there was a certain pattern of non-dipnlc sources that explained 
some puzzling obsersations that he had made of excursions seen in lava How data from Hawaii. 
These observations were that the excursions seemed .ilwavs to bi in the direction such that 
shallower inclinations of the field were produced. Cox suggested tl at this was caused b\ core- 
surface sources drifting past the longitude of H.twaii. which weic preferentially c'lienied and 
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occurred preferentially at certain latitudes. .\t low latitudes the core surface sources would 
be oriented such that the fields they produced would be pointing outward, whereas at 4 . 5 ® 
latitude the core surface sources would produce field pointing downward. Cox ixtinted out that 
tills arrangement would also produce a far-sided effect, which had been noticed by Wilson & 
•\de-Hall (1970), which could, liowevcr, also be explained by an off-centred dipole displaced 
towards the North Poh\ 



Pir.cRE J. .Vormalizcd >urfacc power as a lunciion of decree of harmonic. 

See text for detailed description. 

Unfortunately, the Icelandic data do not support this model. Harrison & Watkins showed 
1979) that the zonal non-dipole field model of Cox should give a slightly near-sided pole at 
Iceland, whereas the data in fact suggest that even here '.here is a far-sided effect (Saemendsson 
ft al. t j8oJ . In addition, the model of Cox predicted that e.xcunions in the Icelandic data 
should be marked by directions that have a higher inclination than normal. search was made 
through all the data collected from Iceland to see if there were preferential directions with 
higher inclinations that were not part of the normal scatter of directions, and no positive effect 
was seen. .Another factor shown by the Icelandic data is that the proportion nf low latitude 
v.g.ps is very large. This was pointed out by Harrison & Watkins (1977), who showed that the 
pattern of v.g.ps for both eastern and western Iceland data do not not fit into the normal Fisher 
distribution of poles. The main reason is that the distribution has too few poles close to the 
mean pole and too many that are a long distance from the mean pole. 

Harrison iic Watkins 1977) showed that a distribution that fitted the data adequately could 
be derived from a Fisherian distribution for most of the poles, with the rest of the poles being 
randomly distributed over the surface of the Earth. 1 he main purpose of this study by Harrison 
& Watkins 1977I was to compare typical continental lava flow results with those obtained 
from lava flows sampled by the Deep Sea Drilling Project, and so the data selection was 
designed to retain d.ita from lava Hows that gave a large amount of within-flow scatter. .\ more 
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rigorous si-Icction of lava Hows was carried out by Dodson ( 1980) wlio found that even with the 
reduced number of low-latitude v.g.ps produced by this more rigorous selection, the Fisher 
distribudon still did not adequately fit the data. Rather, he found that the bipola: distribution 
was a much better fit to the data ^see also Onstott 1980). This distribution also has the ad\ antage 
of not requiring inversion of poles in the Southern Hemisphere before a calculation is made. 
The reason why the bipolar distribution fits the data better is that for the same atigular standard 
deviation the bipolar distribution reaches a ma.\imum (i.e. the mode) at smaller angles than 
does the Fisher distribution, and to make up for this there is a larger proportion of poles at 
large distances from the mean. 
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Kiuure 3. \iigntnrnt of sreondars- rigenveciors for Icelandic and Canary Uland d.sia. plotted relative to site 
longitude. Also shown is a set of secondarv eigenvectors from the statistical tests of Dodson tigSob for 
intermediate latitudes. 


The data selection used by Dodson resulted in about one in five v.g.ps lying at latitudes 
within 50® of the equator. Harrison ( 1980 ) found that l7.7°o of the poles lay further than 40® 
away from the mean pole. In a study of over 1000 lava flow's from northern Iceland, Saemund- 
sson ft al. ( 1980 ) found that of the poles lay within .50® of the equator. Numbers obtained 

from an "xtensive collection of lava flows in the Canary Islands show that at this lower latitude 
the number of v.g.ps found within .50® of the equator is onlv about I in 20. much smaller than 
the number for Iceland. 

There is another remarkable result obtained b\ Dodson (1980' from his analysis of the 
Canary and Icelandic results. When fitting data with the bipolar distribution one obtains an 
eigenvector to a :$ x 3 matrix, which shows the mean direction. The two other eigenvectors are 
in a plane at 90 ® to the principal eigenvector and if the data are not perfectly a.\ially symmetric, 
one of these secondary eigenvectors shows the meridian in which the data tend to be more 
clustered. Dodson found that the directions of these secondarv eigenvector, for the western 
Iceland, eastern Iceland and Canary Islands data sets were in fact verv- close, too close to 
have been caused by chance (figure 3 ). 

The model suggested by Dodson 1980) to explain many of the phettomena is that the non- 
dipole field of the Earth can be approximated by radially oriented dipolo close to the surface 
of the core, as first-order approximations to horizontally flowing curretit loops at the core' 
surface. These radial dipoles are prefercntiallv arranged such that there arc more of them pn 
unit area at high latitudes thatt at low latitiules. T)te preferential l.ttiiuditial arrangement <w 
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Figure 3. Alignment of secondary eigenvectors for Icelandic and Canary 
Island data, plotted relative to site longitude. The asterisk marked 
El Is for the eastern Iceland data set, and the asterisk marked WI,CI 
is for the western Iceland and Canary Island data sets. Also shown is a 
set of secondary eigenvectors from the statistucal tests of Dodson (1980), 
for Intermediate latitudes. 
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such that the prohability uf finding a dipole within 30 ' of the pole was 0.5 compared with a 
probability of o. |3 for a random arrangement. This model was capable of explaining the 
orientation of the minor eigenvectors observed in the Icelandic and Canary Island data. .All 
of the nine intermediate latitude experiments gave intermediate eigenvectors within ± 25° of 
the si e longitude (see figure 3 ). 

Doc son's model is also capable of explaining in a general way the fact that there are more 
low latitude v.g.ps from Iceland than there arc from the Canary Islands. The cere surface < 
dipoles were modelled as normally distributed zero mean dipoles, which meant that there was 
a greater chance of producing a dipole large enough to perturb the dipole field by a large 
amount near Iceland than near the Canary Islands because of the greater probability of finding 
a dipole at the correct location. One thing that Dodson did not calculate was the relative 
strength of the field during times when there were low-latitude v.g.ps. This is an important 
. consideration, in that the data from both cast and west Iceland (Wilson et al. 1972; Harrison 
1980} and another data set from Northern Iceland (Saemundsson et al. 1980) show that the 
average intensity of lava flows giving low-latitude v.g.ps is considerably less than the intensity 
for the other lava flows. .Also, the variation of intensity with v.g.p. latitude for these data sets 
does not look like the data or models presented by Williams & Fuller (i98t) or Hoffman (1979). 



Table 2. Opf-oentre dipole moments 


corrected 


N 

eAi 

St 

v.d. 

r.m.s. 

r.m.t. 

rererence 

8-10 0.38 





0.0077-0.0103 

0.0086-0.0215 

1 


0.24S 

0.01 Itt 

U.08T8 

0.0830 

O.0658 

•y 

8' 

U.2S 

0.0155 

U.0875 

0.0830 

0.0958 

2 

8 

0.26 

U.OI24 

0.0813 

0.0770 

0.0891 

2 

8 

0.28 

0.0032 

O.O620 

0.0581 

0.0673 

2 

20 

0.20AI 

-0.0732 

0.1664 

0.1780 

0.2010 

3 

34 

0.2222 

-0.0478 

0.1805 

0.1 028 

0.2208 

3 

20; 

0.54 

0 

— 

0.0112 

0.0847 

4 

2*3 

' 0.54 

0 

— 

0.0090 

0.0680 

4 


Rel'erencci; 1. Lown (1955); -. Alldredge & Hurwiu (1964): 3. .\lldrrdge & Steams (1969): 4. Dodsuti (19S01. 
4 .Mean dipole moment, normalized b>’ dividing by rj. 


Models of off-centred dipoles 

There have been many models of off-centred dipoles as causes for the non-dipolc field. Some 
of these are listed in table 2 . Lowes (1955) described briefly a model which consisted of eight to 
ten radial dipoles placed at a radius of ■). 3 H (normalized to the Earth's radius), in which the 
dipole moments (normalized by dividing by the cube of the Earth’s radius) ranged from 0.0077 
to 0 . 0193 . The four models by .Alldredge & Hurwitz 1964 were attempts to model various 
forms of the Earth’s magnetic field. They placed the dipoles between <». 24.5 and 0 . 28 . The 
mean dipole moment is quite small because there were both positive and negative dipoles in 
their models. The Itetter measurement of the strength of the dipoles is the r.m.s. value of dipole 
moment. The two models by .Alldredgc & Stearns 1969 used, in our opinion, too mans 
dipoles. The rcfult is '.hat many of them are \er\ close to neighljoun, and positive and negative 
dipoles tend to cancel out. Fhe result is that tne r.m.s. value of the dipole moment for these 
models is considerably higher than the others listed in table 2. Dodson 1980 modelled the 
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non-dipolc field statistically by usini^ 20 dipoles placed at the core acittal boundary. The 
dipoles were randomly chosen from a distribution with zero mean and with a certain standard 
deviation, which is shown in the table. 

All the models except those by Dodson (1980) placed the dipoles well within il;r core, the 
supposition being that by so doing it is possible to make a better approximation to an extended 
horizontal current distribution at the pfe-niantle Ixiundary. We attempt to show that this is a 
reasonable method by calculating the potential produced by the two systems, .\lldredge ( 1980) 
has given formulae for the zonal harmonics g" produced by current loops and also by vertical 
off-centred dipoles placed beneaili the North Pole. We have calculated the variation of magnetic 
potential for current loops Howing at the core mantle boundary as a function of loop radius. 





Fioum 4. Poirniial from current loop*, (a) \'ariation of nomialurd poirnual at Earth’* surface for core surface^ tvA* 

current loopt a* a function of distance from axis. The radius of the loop is ^iven in kilometres, (b) Half 1\ 

width of normalized potential as a function of radius of current loop, (c) Maximum potential of current ' 
loop as a function of radius of current loop. ^ 


These curs es are shown in figure 4 . Note that in the formtilat’-jn of this problem by .Alldredge 
t98o) the centre of the current loop was kept at a constant distance from the centre of the 
Earth. But one consequence of increasing the radius of a current loop flowing at the core mantle 
boundars' is that its centre has to get closer to the centre of the Earth. The 4 ower cun es on thii^ 
figure ^show the magnetic potential as a function of radial distance from the North Pole, for 
current loops of radii between luoo and 2sou km. The cunei are normalized such that the 
maximum potential for the North Pole is unity. The cune isb ov e t hesc]^ls the angulv. half 
width of the potential curves as a function of loop radius. . 4 s would be expected, the wider 
loops give potential functions of greater half width. The cun e^shows the value of the 

potential function at its maximum (North Pole) for unit normalized magnetic moment. Again, 
as would be expected, if the magnetic moment is highly concentrated into a loop of small 
dimension, the maximum value of the potential is large*". 

We have also done similar calculations for vertical off-centred dipoles, using information 
given by Alldredge 1980), Results are presenied in ugure 3 . .\gain, it can l>e seen that as the 
dipole approaches the surface of the Earth increasing vaiue of the half width of the 

potential function at the Earth's surface decreases, and the maximum value of the potential 
increases. We have chosen to model dipoles ,ind core surlace current loops that give the same 
half width of the potential fnnetion by fittimj the two half-width eurM*' In fniure* 4 .ind 3 by 


(UL 


ck\ f^(_ 
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rxpotciitial functions, it is possible to J.ciive an expression for the radius ofilie current loop as 
a function of the depth of the dipole. This expression is 

/C/km = :{ 7 (M - (8) 


Radii calculated for the depths given b\ .Mldredge & Huruitz (1964) range between JoOO 
and 2300 km. We can also use information presented in figures 4 and ■> to determine the dipole 
moment of the current loops. It turns out that »he current loops have a slightly greater dipole 
moment than do the dipoles. This has been indicated in table 2 by showing corrected r.m.s. 
dipole moments, in which tliis conection factor is applied to the r.m.s. values for the dipoles 
themselves. .Although the corrections are only about li",,, this is not necessarily a trivial 
amount, for Cox's (1968) model of field reversal calls for a reversal w hen the axial component 
of the net sum of the off-centred dipoles exceeds th; axial component of the central dipole. 
This probabilitx' is ver\- critically dependent on tlie relative strengths of the central dipole 
and the off-centred dipoles. 






Figcre .3. Same as figure 4. but for vertical dipoles at a distance r, from the centre of the Earth; 

r^ IS the radius of the Earth. 


Dodson (this symposium) emphasizes that the choice of 20 radial dipoles Dodson 19.S0) is 
somewhat arbitrar\-. The important parameter that controls the amount of secular variation 
and the number of low-latitude poles is the total variance of the radial dipoles r,, where 
j. 07 = no7, n is the number of dipole^jand 07 is the variance of an individual dipole. If the r.m.s. 
values for Dodson's dipoles are altered to be correct for eight radial dipoles, the values are 
0.0177 and 0 . 0142 . .A further correction needs to be made to the dipole strengths of Dodson 
1980 Lnd this symposium) because he placed the dipoles at the core surface. If- the dipole 
moment is adjusted for a deeper placement bv using an inverse cube rule, the moments become 
0.0737 and O.Ooftl for dipoles placed at a radial distance of 0.20 of the Earth's radius. .A further 
correction of 15 “0 lo dipole moment of a current loop brings these values up to 0.11S47 

and fi.OtisO, listed in the penultimate column of table 2 . 

This discussion suggests that an appropriate model for the non-dipole sources of core field 
consists of about eight current loops of radius about 24 uo km located at the core surface, and 
having normalized r.m.s. dipole moments lying between tMl 7 and n. 1. 
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C’ONCM'SIONS 

\Vc liavr shovMi that tlir usr of spherical harmonics of low degree to remove the core field is 
not perfect. The field left alter tlie ’core’ field has been removed has sufficient power at long 
wavelengths to make it difficult to imagine crustal sources with enough magnetization to 
produce the desired field satiation. Instead we suggest that off-centred dipoles located within 
the core are a better way to describe tite non-dipole field of ilie Earth. We believe that it is 
possible to remove long-wa\ elength signals more completely by using off-centred dipoles l)oried 
deep within tlie core than Its using spherical ha.inonics. These dipoles represent core surface 
current circulation systems. Data from Iceland and the Canary Islands shoss that viable models 
of the Earth’s non-dipole field can be produced that explain mans of the field variations 
recorded itt extensi* e lava flow piles in these tsvo islands. 

Research was supported by N..-\.S..-\. contract number N.\Sa-2tj20f for work on Migsat 
data. \Ve th~nk Bcb Langel of N..A.S..\. for pros iding us the Gauss coefficients of his degree 2:< 

Magsat field. \Ve have appreciated discussions svith John Southam. Contribution from the ^ 

Roscnstiel School of Marine and .\tmospheric Science, University of Miami. 
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Discusston 

K. M. Creek Department oj Ceophusirs. L nirersitii oj Edtuburgh, L'.k. . Models sucli as that of 
.\lldredge 4 Hurwitz 1904) in w liich liic geomagnetic field is fitted bv several dipoles locatt il \ 

in the outer core plus a central dipole assume the field to be static. That is to say ilic core i' 
assumed to he transparent so iliai fields prodnred b\ snurres 1 w heiher ifie\ be dipoles ttr 1 iirrent 

l‘IS2 fe., 3^, Zso-L^i. 



ORIGINAL PAGE IS 
OF POOR QUALITY 

CORL MAGNI.TK FIl'.I.I) OF IMF. F..\RTH oi.l 

loo|>s ill ill*' r.ir side of the core from ihr observer are transmitted across the core without 
atteiuiatioii. In fact, the i(eom;u>tietic held is not static so the assumption of no attenuation is 
not correct. So my question is; In the models that Dr Harrison has described, has he allowed 
for attenuation of the sii;nal? If not. can lie simijest how this minht be done, beariin; in mind 
ihe fact that the outer core is not solid? 




Reply. We have not allowed for attenuation of the signal. The sources 
are believed to be near the core-mantle boundary because deeper tiue- 
va»y;ng sources would be shielded by the conducting core from having an 
effect at the Earth's surface. This suggests that one possible method 
of treating this problem would be to disregard the effect of any individual 
dipole at locations greater than 90° away from the dipole (l.c. in the 
opposite hemisphere to the dipole). 


